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Abstract 

The second and fourth Fourier harmonics of the azimuthal distribution of particles, V2 and V4, 
have been mesured in Au-Au collisions at the Relativistic Heavy Ion Collider (RHIC). The ratio 
V4KV2) 2 is significantly larger than predicted by hydrodynamics. Effects of partial thermalization 
are estimated on the basis of a transport calculation, and are shown to increase the ratio by a small 
amount. We argue that the large value of V4/(v2) 2 seen experimentally is mostly due to elliptic 
flow fluctuations. However, the standard model of eccentricity fluctuations is unable to explain 
the large magnitude of V4/(v2) 2 in central collisions. 



1. Introduction 

The azimuthal distribution of particles emitted in non central nucleus-nucleus collisions at 
RHIC is a good tool for understanding the bulk properties of the matter created during the colli- 
sions. Near the center of mass rapidity, it can be expanded in the following Fourier series: 
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where <p is the azimuthal angle with respect to the direction of the impact parameter. The large 
magnitude of elliptic flow, V2, suggests that the matter created in Au-Au collisions at RHIC 
behaves like an almost perfect fluid. However, ideal hydrodynamics predicts V4 = |(v2) 2 Q- 
while recent experiments [3, 4] find V4 (vi) 2 . In this talk, I investigate this discrepancy. 



2. Fluctuations in initial conditions 

2.1. Initial eccentricity fluctuations 

Figur^U (left) presents a schematic picture of a non central heavy-ion collision (HIC). The 
overlap area between the colliding nuclei has an almond shape, which generates elliptic flow. 
This shape is not smooth: positions of nucleons within the nucleus fluctuate from one event to 
another, even for a fixed impact parameter. Therefore, the participant eccentricity, epp, which is 
the eccentricity of the ellipse defined by the positions of participating nucleons, also fluctuates. 
Since elliptic flow appears to be driven by the participant eccentricity [5], eccentricity fluctua- 
tions translate into fluctuations of the flow coefficients V2 and V4. 
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Figure 1: (Color online) Left: Picture of the two frames used for defining the initial eccentricity (from [5s]). The x 
axis defines the reaction plane while the x' axis is the minor axis of the ellipse defined by the participants (grey dots). 
Right: Centrality dependence of v 4 /(v 2 ) 2 : data from STAR |6] and PHENIX |7]; eiTor bars on STAR data points are our 
estimates of nonflow errors 1 1]. Lines are predictions from ideal hydro with or without fluctuations. 

2.2. Impact of flow fluctuations on v 4 /(v2) 2 

There is no direct way of measuring vo and V4. Analysis methods rely on multiparticle cor- 
relations. Experimentally, V2 can be extracted from the 2-particle correlation and V4 from the 
3-particle correlation using (cos(20i - lipi)) = ((V2) 2 ) and (cos(4</>i - 2<pi - 2^3)) = (v4(v2) 2 ), 
where angular brackets denote an average value within a centrality class. Thus any experimen- 
tal measure of V4/(v2) 2 , obtained using these methods, is rather a measure of (V4(v2) 2 )/((v2) 2 ) 2 . 
Inserting the prediction from hydrodynamics v 4 = j(v2) 2 , we obtain 



We assume that v% scales like epp, whose fluctuations can be estimated using a Monte-Carlo 
Glauber model [5]. The resulting prediction for V4/(v2) 2 is displayed in figure Q] (right). Fluc- 
tuations clearly explain most of the difference between hydro and data. However, experimental 
data are still slightly higher than our prediction from fluctuations. We argue that for peripheral 
to midcentral collisions, the small residual difference may be understood in terms of deviations 
from local equilibrium. 

3. Partial thermalization effects 

So far, we have assumed that ideal hydrodynamics correctly describes the expansion of mat- 
ter created in a HIC. But ideal hydro assumes that the system remains in local thermal equilib- 
rium (regime where the average number of collisions per particle n co u is large) throughout the 
evolution. In a previous work 18|] we have shown that, in order to reproduce the centrality depen- 
dence of elliptic flow, the deviation from local thermal equilibrium must be taken into account 
(ftcoii K 3 - 5 would be a typical value for Au-Au collisions at top RHIC energy). 

In the limit of small n ca u, one expects both V2 and V4 to scale like n co u, so that V4/(v2) 2 scales 
like 1 /n co u: we thus expect that the farther the system from equilibrium, the larger v 4 /(v2) 2 H- 
In order to have a more quantitative estimate of this effect, we use a 2+1 dimensional solution 
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Figure 2: (Color online) Left: Variation of V4/(va) with the Knudsen number. The point at K = are obtained using an 
independant ideal hydro calculation. Right: Same plot as figure[T](right), with one additional curve showing the effect of 
the deviation from local equilibrium. 



of the relativistic Boltzmann equation to study systems with arbitrary n co u. We use the Knudsen 
number [9], K oc l/n co u, as a measure of the degree of thermalization of the system. Figure 2 
(left) displays the variation of V4/(v2) 2 with K (see [1] for details). Extrapolation to the hydro- 
dynamic limit K — yields the value 0.52, quite close to the expected I. For nonzero values of 
K, V4/(v2) 2 slightly increases. The effect of this increase on the centrality dependence is shown 
in figure 2 (right). The values of K are borrowed from a previous study [8]. When both fluc- 
tuations and partial thermalization are taken into account, our calculation slightly overshoots 
data for midcentral and peripheral collisions, but the overall agreement is good. We do not yet 
understand the large value of V/t/(v2) 2 for central collisions. 

4. Conclusion 

We conclude that: 1) v\ is mainly induced by vi\ 2) the deviation from local equilibrium has 
a small effect on V4/(v2) 2 ; 3) eccentricity fluctuations explain the observed values of V4/(v2) 2 , 
except for the most central collisions which require further investigation. 
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